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We present two methods for determining the absolute detection efficiency of photon-counting de- 
tectors directly from their singles rates under illumination from a nonclassical light source. One 
method is based on a continuous variable analogue to coincidence counting in discrete photon ex- 
periments, but does not actually rely on high detector time resolutions. The second method is 
- - - based on difference detection which is a typical detection scheme in continuous variable quantum 

, optics experiments. Since no coincidence detection is required with either method, they are useful 

3r ■ for detection efficiency measurements of photo detectors with detector time resolutions far too low 

' to resolve coincidence events. 

CN ■ 

^ ■ I. INTRODUCTION 

<■ 

' High detection efliciency is crucial in many contexts within quantum optics. In particular, recent linear optics 
quantum computation methods have been shown to require detection efficiencies higher than 90% to be scalable P, 
' y> Hi S ' Another well known example comes from Bell inequality experiments where a detection efficiency of at least 
J> i 2/3 is required to close the detection loophole 0. Thus higher detector efficiencies do not only have the obvious 
advantage of providing more signal, but are also of great relevance for both current quantum information science 
and fundamental quantum physics. Therefore building new high efficiency detectors is a crucial area of research. In 
this article we present methods to evaluate their absolute detection efficiency. The typical way to measure detection 

■ efficiency is to use a calibrated reference detector, and to compare its response to an incident light beam of constant 
\^ , intensity with the response of the detector under test. It has been shown that quantum mechanics itself provides a 

' way to measure the absolute detection efficiency of light detectors without the need for a reference detector. By using 
non-classical photon statistics such as those produced in spontaneous parametric down conversion 6] both theory 
^ ^ and experiment show that the absolute detection efficiency can be determined by the ratio of the coincidence rate 
I \ to the singles rate 013, 0- However, this method requires the detectors to have a high enough time resolution for 
coincidence counting. We provide work-around schemes that overcome this limitation. Specifically, we have developed 
two related detector efficiency calibration methods that use quantum correlated light but do not rely on high time 
^ . resolution coincidence counting |0 . This might be of particular interest for novel detectors and prototypes which are 
^ '' not yet capable of coincidence counting, such as electron multiplying CCD cameras. 

^ . The first method is based on a continuous variable analogue of coincidence counting in discrete photon experiments. 
Recall that coincidence counting is essentially bitwise multiplication (i.e., an AND gate) of single counts within a small 
. , coincidence window. In the continuous variable limit, this is achieved by measuring the mean product of the detected 

■ photon numbers in two beams generated by spontaneous parametric down conversion. Since down conversion photons 
are emitted in pairs, the normalized mean product of the photon numbers measured in the two beams during a specific 
time is maximal in the case of perfect detection efficiencies when all photons are detected. On the other hand, in the 
limit of small detection efficiencies, the normalized mean product of the detected photon numbers has its minimum. 
For the general case an expression can be derived, that allows calculation of detection efficiencies from the measured 
mean product of the detected photon numbers in the two down conversion beams. 

The second method uses the variance of the detected photon number differences in two beams generated by spon- 
taneous parametric down conversion. Such difference detection is a typical detection scheme in continuous variable 
quantum optics experiments 11, 12]. Since the photons in the two beams are perfectly correlated, the difference in the 
photon numbers measured in the two beams is zero in the case of perfect detection efficiencies when all photons are 
detected. Uncorrelated loss in these two beams diminishes those perfect correlations. Therefore in the limit of small 
detection efficiencies, the normalized variance of the differences in the detected photon numbers has its maximum. 
For the general case again an expression can be derived, that allows to calculate the detection efficiencies from the 
measured variance of the detected photon number differences in the two down conversion beams. 

The article is structured as follows: First we derive these relations mentioned above disregarding background. 
Then we generalize these results to include two different background levels in the two down conversion beams, since 
background light is a very significant contribution in single photon counting experiments. 
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II. THEORY 
A. Product Detection Method 

We begin by deriving a relationship between the mean product of the singles rates and the detection efficiencies in 
each down conversion beam, rji and 772, in the absence of background. We assume that parametric down conversion 
emits light beams described by a general distribution G'Ar(fc) of the number of photon pairs k with mean value 
{k)G = N and the second moment {k'^)G E3- The probability of detecting / out of k photons in each of the two down 
conversion beams is given by the binomial distribution Bk^rnil) = (/) vi (1 ~ f^i)*^ '1 with mean value {1)b = rjik and 
the second moment {P)b = Vi^ ~ vfk + 'nfk'^ (* = 1, 2). Thus the mean product of the detected photon numbers in 
the two down conversion beams is given by 



00 k k 

k=o 1=0 m=a 



The corresponding result for 772 is 



= mv2{k')G = ^i{i')-{i) + m{i)), (1) 

where I and m are the detected photon numbers in each down conversion beam, (l) — rii{k)G and {l'^) — rii{k)Q — 

nl{k)G + vl{k^)G. 

Together with the expression ~ the detection efficiency 771 follows as 

^2 = ^ - + 1 . 3 

Note that neither formula depends on the coincidence rate. However, the quantum statistics of the light enters the 
expressions in the mean product of the singles rates. 

In a former method 7| the absolute detection efficiency is determined from the ratio of the mean coincidence rate 
(c) to the mean singles rate (Z) or (m), 771 — and 772 = Together with or J^l one obtains for 77 = 771 =7^2, 
and hence (s) = (l) = {m) and (s^) = (P) = (771^), 

(,s)-(c) = (,s^)-(M^fc;^. (4) 



This simple expression relates the difference between the mean singles rate and the mean coincidence rate to the 
variance of the detected photon number differences in the two down conversion beams. This motivates our second 
approach for determining the detection efficiencies from the variance of the detected photon number differences in the 
two down conversion beams described in section HTH 

B. Product Detection Method: General Approach Including Background 

We now extend this theory to cover more realistic experimental conditions and correct the measurements for possibly 
different backgrounds in the two detectors. The averaged quantities contained in (|2J) and Q have to be extracted 
from experimentally accessible quantities which include background. We do so by splitting up the measured photon 
numbers (subscript M) into the photon numbers corresponding to the signal (no subscript) and into photon numbers 
corresponding to the background (subscript B), where I — Im — Ib and m = vtlm — ms- The background can be 
estimated experimentally from a separate configuration. We get 

{1} ^ ihi) - ih) , (5) 
(777.) = (777m) - {rriB) , (6) 
{I') = (Ih) - Hi) - 2{Im){Ib) + 2{Ib)\ (7) 

(777^) = (777^^) - {ml) - 2 (777Af)(777s) + 2 (777s)' , (8) 
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and 

(Im) = (iMmM) ~ {hi) {ins) ~ {lB){mM) + {lB){'mB) ■ (9) 

Here we used the statistical independence of / and Ib and m and to_b, respectively. By inserting these expressions 
into (0) and Q the detection efficiencies can be determined from the data directly measurable in an experiment: 

^ ihimM) - {lM){mB) - {lB){mM) + {lB){mB) _ {Ijj) - {1%) - 2 (Im){Ib) + 2 {hf ^ ^ ^^^^ 
{itim) - {tub) {Im) - {h) 

and 

{himM) - {lM){mB) - {lB){mM) + {lB){mB) (m^j) - (m|) - 2 (mM)(mB) + 2 (ms)^ 

= Ti — ^ — Tr~\ 1 — ^ — '/ — ^ ^ ■ 

{I.M) - {Ib) {mn) - {tub) 

C. Difference Detection Method 

Our second approach for determining the absolute detection efficiencies in each down conversion beam, jyi and 772, 
relies on measuring the variance of the differences in the singles rates {{I — mf'). This method is closely related to 
the approach described in section Hi Al since {{I — raf') — {P) + {m?) — 2 {Im), and both ((/ — to)^) and {Im) depend 
on the degree of second-order coherence which is affected by uncorrelated loss. However, the two methods may be 
useful under different circumstances, especially since difference detection is a typical detection scheme in continuous 
variable quantum optics experiments. 

The variance of the detected photon number differences in the two down conversion beams is given by 

OG k k 

{{l~mf) = Y,Y.ll G^(fc) ^'^.'u (0 Bk,r,, (m) (/ - mf 

k=0 1=0 m=0 

= {I) + (m) - {m{l) + ,?2(m)) + , (12) 

Vi 

where I and m are the detected photon numbers in each down conversion beam. 

_W_ = 22i 

(m) Jj2 ' 



Together with the expression 4^ = the detection efficiency 771 follows as 



Correspondingly, 



3{m)-% + {P)il-<^)^-{il-m)^) 
2 {m} 



" 2(0 ■ ^ ' 



For 77 = 771 = 772, and hence (s) = (/} = {m), 11211 becomes 

((^-7r^)^)=2(l-7;)(.). (15) 

For perfect detection efficiencies (77 = 1), ((/ — 777,)^) = 0. Since the down-converted photons are always created in 
pairs, the difference in the detected photon number, and hence the variance of that quantity, is exactly zero if all of 
the pairs are detected. In the limit of very small detection efficiencies (77 -C 1), {{I — m)"^) k, 2(s), which corresponds 
to the variance of two independent Poissonian light beams of equal intensities. The non-Poissonian contributions in 
G^ik) cancel out, showing that quantum statistics strongly depends on the detection efficiency of the detectors. 
Expression \ib\ can be rewritten as 

A similar expression was stated without explicit derivation previously 0, 0| . 

With (s) = r]N, the normalized expression (|15|l reads ^ = — 1), which verifies that the normalized 

variance of the detected photon number differences in the two down conversion beams diverges for 7; — > and goes 
to zero for 7; — > 1. Using equation I^J, an analogous treatment can be performed for the normalized mean product of 
the detected photon numbers in the two beams. 
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D. Difference Detection Method: General Approach Including Background 



As in section rilBI wc extend this theory to cover the more reahstic experimental conditions including background. 
The quantities in H13|l and 114|) have to be extracted from quantities that are directly accessible to measurement. We 
get ©-(IHl) and 

{{1 ~ mf) = {{hi - ruMf) + 2 {{Ib) - (ms))' - 2 {{Im) - {mM)){{lB) - (ms)) + 2 {lB){mB) - {1%) - {m%) (17) 

for the background corrected difference term. 

By inserting (jSJ-® and H17|l into and 114|) the detection efhciencies can be determined from the data directly 
measurable in an experiment, 

^ H{mM)-{mB))- ''^"lrJ)~Xf ^ {llt)-{ll)~2{lM){lB)+2{lB? A _ {vn,) - {^bW 
VI ^ ^^ + o./„. ^ ^^ {Im) - {h) J 

(18) 



2{{mM)-{mB)) 2{{mM)-{mB)) \ {Im) - {Ib) J 

{{hi - niM)^) + 2 {{Ib) - (mg))^ 2 {{hi) - (mM))((^B) - (mg)) ^ 2 {lB){mB) - {1%) - (m|} 

2 ((tom) - (™b)) 2 ((mAf) - (™b)) 2 ((mM) - (ms)) 

and 



3((/Af)-(/B))- iSJg^) - {ml) - 2 (mM)(mB) + 2{mB)^ ( _ {hi) - {h) 

2{{hi) ~{Ib)) 2{{hi)~{lB)) V (™M)-(mB), 

{{Im - mu)^) + 2{{Ib) - {tiib)? 2 {{Im) - {mM)){{lB) - (ms)) , 2 (/s)(mB) - (/|) - (m|) 



2((?Af) - (^b)) 2((Zm)-(;b)) 2((?m)-(Zb)) 



(19) 



E. Error Estimates 



Finally, we want to derive the statistical errors of the detection efficiencies (without background), i.e., {A) for the 
product detection, eq. {B) for the difference detection, eq. l(T^ . and (C) for the coincidence method {7^, described 
above eq. 

.1^^ = 14-^ + 1, (20) 
(m) (/} 

(B) 3(m)-<^ + (0(l-illf)^-((^-m)^) 

= 2M ' ^ ^ 

77^^) - (22) 
" (m) ' ^^^^ 

In each case rji is a function of several mean values, i.e., iji = r/i{{u), {v), ...) where u,v, ... G {l,m,l^,lm, {I — m)^,c}. 
The sample variance of 771 is defined as 

^ (^i) = ^u>(^aRj +''<^na(^j ^^^^-^^''^aM + - ^^^^ 

where 

^(u> = M ' ^^^^ 

_ {uv) ~ {u){v) 

'^{u){v) = (25) 

are the variances and covariances of the sample means with sample size M. It has to be stressed that in experiments 
the time interval chosen for accumulating the individual measurements needs to be much larger than the resolving 
time of the detector under test. The mean values are given by 



fc fc 

fc=0 1=0 m=0 



) = E E E (^) ^^^''1 (') ^^'V. (m) X . (26) 
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FIG. 1: Sample variances cr^(ryj'*-') (product method, solid line) and cr^(»7i^^) (difference method, dashed line) as a function of 
r]i for a Poissonian distribution in the limit TV 3> 1. For the sake of generality the sample size is = 1, where all variances scale 
with 1/M. (a) The second efficiency is constant: r]2 = 0.1. The variance cr'^iVi^^) diverges for 771 0. (b) Equal efficiencies: 
J72 — rji. 



Due to the perfect correlations the first and second moment of the coincidences can be computed by applying the 
binomial distribution to one down conversion arm twice: 



k I 



(^') = Gjv(fc) Bk,,, (l) Bk,,, (m) mP , (27) 

fc=o ;=o m=0 

where p = 1,2. 

For a Poissonian down conversion distribution Gjv (k) — N'' exp(— A^) /kl, with N the expected mean photon number 
in one sample measurement, the efficiency sample variances in the three cases are 

2f (A). Vi{l + N- f]i) +Nr]2[2 + m iVi - 4)] 

"('^i ) = Wmf, ' 

2 M iV 77J ?72 

In the limit N ^ 1 the variance for the coincidence method vanishes: cf'^{'n{^^) — > 0. For the product and difference 
detection the variances approach constant values which depend on the efficiencies 771 and 772 as well as the sample size 
M . Figure 1 shows these two variances as a function of 771 for M = 1 where all sample variances scale inversely with 
the sample size. In general, the difference method is more accurate than the product method, except in the case of 
fixed 772 and vanishing 771. In the special case of equal efficiencies 771 — 772 the expression for the difference method 
simplifies tremendously: 

(31) 

For a thermal distribution Gjv(fc) cx exp(— fc/A'^) the first two variances are more cumbersome and we do not write 
them here. In the limit of increasing N the variance for the coincidence method a'^{r][^^) vanishes again, while the 
variances for the product and difference detection method linearly diverge in the limit TV 3> 1, still also scaling with 
1/M. Only in the special case of equal efficiencies 771 = 772 the variance for the difference method a^{ri[^^) becomes 
independent of TV. Therefore, if equal detectors are used, the difference method is more favorable. In this special case 
and for TV ^ 1 we have the simple expression 

(32) 



Hence, in this case, the variance for the thermal distribution has the same form as expression H31|l for the Poissonian 
distribution. 
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III. CONCLUSIONS AND OUTLOOK 



We have presented two methods for determining the absolute detection efficiency of photo detectors. Since they are 
apphcable to detectors with low time resolution, they overcome the limitations typical for absolute detection efficiency 
measurements. The first is based on measuring the mean product of the detected singles rates in two beams generated 
by spontaneous parametric down conversion. The second method uses the variance measurements of the differences in 
the detected singles rates in the two down conversion beams. The two methods correspond to the different detection 
methods typically used in either the discrete photon or continuous variable communities, respectively. Both procedures 
could be used for measuring the absolute detection efficiency of photo detectors that do not provide the appropriate 
time resolution for coincidence counting. 
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